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Abstract. We present statistics (S-statistics) based only on random variable 
(not random value) with a mean squared error of mean estimation as a concept 
of error. 
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Let us consider (e.g. for j = n + 1) variance of the difference Rj of two random 
variables Vj and Vj, where E{Vj} ~ E{Vj} ~ m, in terms of covariance 

D^V,-V,} = Cov{{V,-Vj){V,-V,)} 

= Cov{V,V,} - Cov{V,V,} - Cov{V,V,} + Cov{V,V,} 
Cov{V,V,} ~ 2Cov{V,V,} + Cov{V,V,} 



2 



T. SUSLO 



introducing the estimation statistics Vj — ^jVi — ^jVi 

D^{Rj} = Cov{VjVj} -2Cov{VjVj} + Cov{VjVj} 

= Var{V,} - 2Cov{J:^ lo^V^V^} + Cov{{Y., ^iVi){Y.. ^)Vi)} 
- 2 Ei '^]Cov{V,V, } + E; Lo]u:]Cov{ViVi} 

in terms of correlation function p 

D'^{Rj] = CT^ + 2(72^^^^. - a'^ujipiifj) 
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The unbiasedness constraint (I condition) 

E{Rj} = E{Vj - Vj} = E{Vj} - E{Vj} = 

is equivalent to 



The minimization constraint 



where 



produces n equations in the n + 1 unknowns: kriging weights uij and a Lagrange 
parameter fxj (II condition) 
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multiplied by 
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and substituted into D'^{Rj} 

D'{R,} = E{[V,-V,n-E'{V,-V,} 

V ' 

Em^m)-iV,-m)n 
= E{[Vj - m]2} - 2{E{V^Vj} - m^) + £;{[t>j - mf) 
- 20-2 Iwipy I + cr2 |w]pi,a;i | 

give the minimized variance of the field Vj under estimation 

E{{{Vi - m) - - m)f} = a'{l ± (a;}^,,- + /x,)) 
and these two conditions produce n + 1 equations in the n + 1 unknowns 
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2. Asymptotic properties of origin 



Since 



then 
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and (since) 
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the minimized variance of the field Vj under estimation 

- m) - {V^ - m)f} = a\l ± {uj]pi^ + 
has known asymptotic property 

lim E{[Vj - V^f} = lim E{[Vj - uj)vif} = E{[Vj - mf} = 

3. The field under estimation 
Let us consider the field Vj under estimation 



Vj = LO'jV , 



where for auto-estimation holds 



v = 5iv = v 
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with minimized variance of the estimation statistics 

E{[Vi-mf} = _Ccyv{{ujiVi){ujiVi)} 



where for auto-estimation holds 

E{[Vi-mf} = E{[Vi-mf} = a^ 

that means outcoming of input value is unknown for mathematical model, with 
minimized variance of the field Vj under estimation 

E{[(V,-m)-(Vj-m)f] = ^^(l ± (^p,,- + M,)) 
where for auto-cstimation holds 

E{[{Vi - m) - (V- - m)f} = E{[V, - mf} - 2{E{ViVi} - m^) + E{[Vi - mf} = 

^ V ' ^ V ' ^ V ' 

(72 2(72 (72 

that means variance of the field is equal to variance of the (auto-)estimation sta- 
tistics (not that auto-estimation matches observation). 



For j ^ 00 



4. Generalized least squares solution 
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and a disjunction of the minimized variance of the field Vj under estimation 
Em - m?} - {E{VjVj} - m') + E{Vj[Vj - Vj]} if pijcj^ + ^ . = ^ + ^ . = 

^ V ' ^ V ' 

which fulfills its asymptotic property the kriging system 
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equivalent to 
and 

F'uj = 1 

where: F = (1, . . . , 1)', w = (w^, . . . A = A' = pa; i = 1, . . . ,n, has the 

least squares solution 

1 

and 
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with a mean squared error of mean estimation 

E{[V, - mf} = T^^2^ ■ 
5. Approach to simple statistics (conclusion) 
For white noise 

E{[{Vj-m)-{Vj-m)]^} = E{[Vj - m]^} + E{[V, - m]^} 
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Um E{[Vj - m]2} = 
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then 
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and 
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Remark. Precession of arithmetic mean can not be identical to cause a straight 
hne fitted to high-noised data by ordinary least squares estimator can not have the 
slope identical to 0. For this reason the estimator of an unknown constant variance 

1 sum of deviation squares 

n 

in fact is the lower bound for precession of the minimized variance of the field under 
estimation 

E{[{Vj-m)-{Vj-m)f} = a' (l + 1) 

to increase 'a bit' the lower bound for n ^ cxd we can effect on weight and reduce 
total counts n by 1 because 1 is the closest positive integer number to so it is 
easy to find the closest weight such that 

1 1 

> 



n — 1 n — 
then the so-called unbiased variance 

1 



sum of deviation squares 



n — 1 

in fact is the simplest estimator of minimized variance of the field under estimation. 

References 

[1] E. H. Isaaks and R. M. Srivastava, An Introduction to Applied Geostatistics, New York: 

Oxford Univ. Press (1989) 279-307. 
[2] http://www.compsoc.man.ac.uk/~lucky/Democritus/Theory/corfunl.html 



